Quantum state tomography with quantum shotnoise 
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We propose a scheme for a complete reconstruction of one- and two-particle orbital quantum 
states in mesoscopic conductors. The conductor in the transport state continuously emits orbital 
quantum states. The orbital states are manipulated by electronic beamsplitters and detected by 
measurements of average currents and zero frequency current shotnoise correlators. We show how, 
by a suitable complete set of measurements, the elements of the density matrices of the one- and 
two-particle states can be directly expressed in terms of the currents and current correlators. 
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According to the standard interpretation of quantum 
mechanics the wavefunction, or more generally the den- 
sity matrix, determines the probabilities for the possible 
outcomes of any measurement on the quantum state. To 
completely characterize the wavefunction of the state is 
therefore of fundamental interest 0. It is however im- 
possible to infer anything about an unknown state from a 
single measurement, a complete characterization requires 
an ensemble of identically prepared states and the mea- 
surement of a complete set of observables on the state 
Q ■ A reconstruction of the quantum state wavefunction 
via such a series of measurements is known as Quantum 
State Tomography (QST) 0. 

Initially, QST was performed experimentally on the 
discrete angular momentum state of an electron in an hy- 
drogen atom 0. During the last decade QST has been 
performed on e.g. the quantum state of squeezed light 0, 
the vibrational state of a molecule || , the motional state 
of trapped ions Q and of atomic wavepackets [j|. Re- 
cently there has been an interest in QST of two-particle 
states in the context of quantum information process- 
ing. The entanglement of a quantum state, a potential 
resource for quantum information processing, is charac- 
terized by the density matrix of the state. The quantum 
state of polarization entangled pairs of photons has been 
reconstructed using QST 

To date, no QST has been performed on quantum 
states in solid state systems. Very recently a theoreti- 
cal scheme was developed for solid state two-levels 
systems, qubits, appropriate for e.g. the macroscopic su- 
perposition state in superconducting qubits and the spin- 
state of electrons in quantum dots. The set of measure- 
ments necessary to reconstruct the state involves con- 
trolled rotations and detection of the individual qubits. 
For coupled qubits, where entanglement between the 
qubits is of interest, such measurements are highly in- 
volved and have not been demonstrated. 

In this paper we take a different approach and present 
a scheme for QST of discrete single and two-particle or- 
bital quantum states in mesoscopic conductors. The or- 
bital quantum states [IlIIT3 . ll3j | are continuously emitted 



from the conductor during transport, making a long time 
measurement equivalent to an average over an ensemble 
of states. The orbital states can be manipulated by elec- 
tronic beamsplitters, experimentally available [lil 
and detected by measurements of average currents and 
zero frequency current correlators, shotnoise [Til 
This scheme, with all components experimentally real- 
izable in e.g. Quantum Hall systems |l5l Il8|. allows for 
a complete characterization of the quasiparticle quantum 
state in mesoscopic conductors. 

The key question for any QST is: what quantum states 
with interesting properties can be investigated with ac- 
cessible experimental technics? In mesoscopic conduc- 
tors, one typically measures electrical currents and cur- 
rent correlators. In several recent works [ill [Til fl3l IT9L 
lifij . it has been shown theoretically that quantum cor- 
relations, entanglement, between two spatially separated 
particles can be investigated via current correlation mea- 
surements. In particular, in Refs. [ill ll-'il| it was shown 
how entangled orbital quasiparticle states could be gen- 
erated, manipulated with experimentally available elec- 
tronic [T3 . Hil l beamsplitters and detected via current 
correlation measurements. Orbital one- and two-particle 
states investigated via current and current correlations 
are thus natural candidates for QST in mesoscopic con- 
ductors. 
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FIG. 1: Schematic of the setup. A mesoscopic conductor 
acting as a source S (red box) is connected via four leads, 
Al, A2, Bl and B2, to regions A and B (dashed boxes), each 
containing two reservoirs, + and — , a beamsplitter (green 
cross) and a sidegate (blue open box) which induces a phase- 
shift (f>A or <f>B- Two particles (red dots) emitted from the 
source are shown. 
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A generic setup for such orbital QST is shown in Fig. 
^ A mesoscopic conductor S acts as a source for or- 
bital quantum states. The source is connected via four 
single mode leads A1,A2,B1 and B2 to two regions, A 
and B, where the emitted state is manipulated and de- 
tected. The mesoscopic source has one or more reser- 
voirs biased at eV and an arbitrary number of reservoirs 
kept at ground. We note that two-particles effects are 
only present for two or more reservoirs biased The 
temperature is taken to be zero. It is assumed that the 
scattering in the conductor is elastic, however arbitrary 
debasing inside the conductor can be accounted for. 

The regions A and B each contain an electronic beam- 
splitter 0, 0] and an electrostatic sidegate (see e.g. 

to induce a phaseshift, (f>A or (f>B, by modifying the 
length of the lead. The beamsplitters, taken to be reflec- 
tionless, are further connected to two grounded reservoirs 
+ and — where the current is measured. The combined 
beamsplitter-sidegate structure can be characterized by 
a scattering matrix, for e.g. A given by 

v /Xj e i(¥>A3-<£A) \ 

r W7pi(<PAl+<PA3-<pA) I ' V > 
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The transmission probability Ta = 1 — Ra can be con- 
trolled via electrostatic gating [Til Hil Hoj . The phases 
<PAii i = 1, 2, 3 picked up when scattering at the beam- 
splitter are however assumed to be uncontrollable but 
fixed during the measurement. 

The quantum state emitted by the mesoscopic source 
is in the general case a manybody state, it is a linear 
superposition of states with different number of particles 
[21|. However, one- and two-particle observablcs such 
as current and noise are only sensitive to the one- and 
two-particle properties of the state. These properties are 
quantified by the reduced density matrix, which thus is 
the object of interest. Only in some special ca ses, typ- 
ically in conductors in the tunneling limit [111 ll2L , 
are the emitted states true one or two-particle states. In 
the presence of dephasing, the emitted state is mixed. 
Moreover, even an emitted pure manybody state gener- 
ally gives rise to a mixed reduced one- or two particle 
state. It is therefore appropriate to discuss the state in 
terms of density matrices. 

To simplify the discussion we consider a spin-polarized 
system with scattering amplitudes independent on energy 
on the scale of the applied bias eV , i.e. the linear voltage 
regime. The emitted state then has only orbital degrees 
of freedom. We first consider the single-particle orbital 
state emitted e.g. towards A (same considerations hold 
for B). Introducing operators b An creating electrons in 
lead An, with n = 1,2, propagating out from the source, 
the 2x2 density matrix (not normalized) is by definition 
given by 



2 

PA = ^ ' Pmr, 
n,m— 1 



b An \0)(Q\b Am = 



Pn 



P22 



(2) 



where we work in the basis {|1}ai |2}a}, with 6^ n |0) = 
\n)A, formed by the lead indices (see Fig. QJ. The ma- 
trix elements p nm = (b Am bA n )- The Hermitian density 
matrix, pa = p\, has four independent parameters and 
can be written as follows 



PA 



i=0 



Co + c 3 ci - ic 2 
ci + ic 2 c - c 3 



(3) 



where {fTi} = [l,a x ,a y ,a z ]. A normalized density matrix 
is obtained by dividing all elements by cq. In the same 
way, the two-particle density matrix is given by 



PAB= J2 Pnmb A n^Bk\0)(0\b Bl bAr. 



(4) 
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( b Am b Bl b BkbAn)- The 



with the matrix elements p„ m 
two-particle density matrix has 16 independent paramc 
ters and can be written 
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i,j=0 



with £g> the direct product. Expressing the real coeffi- 
cients Ci and Cij in terms of outcomes of ensemble aver- 
aged measurements thus gives a complete reconstruction 
of the emitted state. The accessible measurements are 
average current and zero frequency current correlations. 
Importantly, in the transport state the source continu- 
ously emits quantum states. As a consequence, the long 
time measurements automatically provide an ensemble 
average measurement. At A the average currents at con- 
tacts a = ± are 
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The zero frequency correlator between currents fluctu- 
ations AI in reservoirs Act and B/3, given by S A g - 



(1/2) / dUAI Aa (0)AI Bp (t) 
written [lq j 



AI Bp (t)AI Aa {0)) can be 
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with 



b [ Ba bBa- The operators bAa and bsp in the 



reservoirs at A and B are related to operators bAn and 
bsk in the leads An, Bk, with n, k = 1, 2, (see Fig. ^) via 
the scattering matrix Sa of the beamsplitters at A as 



b A - 



S A 



&A1 
bA2 



(8) 



and similarly at B. 

We start with the reconstruction of the one-particle 
state at A, accessible via the average current (the same 
procedure holds for the state at B). Here a formal ap- 
proach is taken which directly can be extended to the 
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investigation of the two-particle state. We note that the 
reconstruction approach is similar to QST schemes for 
qubits in other systems, see e.g. Refs. 0,H3|. There 
are however a number of important special features for 
mesoscopic systems, making a detailed investigation im- 
portant. It is desirable to minimize both the type and 
number of experiments having to be carried out. As is 
clear from the following, for a complete reconstruction it 
is sufficient to consider only measurements of currents in 
one reservoir in A. Here we consider the current at A+. 
Using the relation between operators, Eq. ©, and Eq. 
ijBJ, we have 



I+/(e 2 V/h) = (n+}=tr(p A A) 



(9) 



with the matrix 



A = 



Ra 

^T A RAe t{,pA+VA) 



^/TlRAe-^+f^ 



(10) 



The phase (fA = 4>A2 — 4>A3 contains all the information 
about uncontrollable phases of the beamsplitter. From 
Eqs. © and i|10|> it is clear that the phase ip A can be 
included in pa by a change of local basis pa — > UaPaUa 
with XJ a = diag[exp(— exp(i^/2)]. Below we 
consider the reconstruction of pa, parameterized by co- 
efficients Ci [see Eq. (0], thus working with A(ip A = 0) 
in Eq. I|1(J|I . This yields pa up to an unknown local basis 
rotation. 

Importantly, only four settings of the beamsplitter are 
needed, both for the current and the current correlators, 
for a complete state reconstruction. The settings / to 
IV considered here are listed in the table in Fig. [21 By 
constructing suitable linear combinations ja{j) of the ob- 
servables at the different settings, in the {|1)a, |2)a} ba- 



Setting 


Parameters 


1 


T A =0, (j> A arb. 


II 


T A =1,<|> A arb. 


III 


T A =1/2, 4> A =0 


IV 


T A =1/2, (|> A =7t/2 




FIG. 2: Left: Table with scattering parameters at A for the 
four different settings. Right: Schematic of an elementary 
mesoscopic source S (red box, see Fig. 0. Two reservoirs 
biased at eV, 1 and 2, and two grounded reservoirs, 3 and 4 
are connected via beamsplitters to the four leads going out 
towards A and B. Scattering between upper and lower leads, 
e.g. Al and A2, is not possible. 



j A (0) = A(I)+A(II) = 1, 

j A (l) = 2A(III) - [A{I) + A(II)} = a x , 

j A (2) = 2A(IV) - [A(I)+A(II)} = <r y , 

j A (3) = A(I)-A(II) = a z , (11) 

we obtain a complete set Q of measurements, since the 
Pauli matrices <jj obey the relation tr(<7j0j) 
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Here A(I) is the matrix A in Eq. (|10|l for the set- 
ting I etc. From Eqs. © and (|llfl and the relation 
(j A {j)} = tr(pA<yj) = Cj we then directly obtain the co- 
efficients Cj, parametrizing pa in Eq. Q 
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in terms of the measured currents for the different set- 
tings, taking the index {k} = [0, 1,2,3] = [I, II, II, IV] . 
This completes the one-particle state reconstruction. 

We then turn to the two-particle state. In Eq. J7J), 
the quantity that is directly linked to the density matrix 
elements p^ m is the reducible correlator (n A n^). This 
correlator is directly obtained from the measured noise 
and the average currents. In analogy to the current, it is 
sufficient to consider correlations between currents in one 
terminal in A and one in B. Considering here A+ and 
B+, one obtains from Eq. (7J and (JSJ the dimensionlcss 
correlator 
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1 A 1 B 



2e*V/h ' p V/h) 2-< n te)=«0>ABA9B) (13) 

where the matrix B is given from A in Eq. I|10|l by chang- 
ing indices A — > B in the scattering amplitudes. Similar 
to the one-particle state, we note from Eq. I|13f) that both 
phases ipA and tps can be included in pab by indepen- 
dent local rotations pab — * {Ua ® Ub)p~ab{Ua ® Ub)\ 
with Ub = diag[exp(— i(^s/2),exp(i<^s/2)]. Below we 
thus consider the reconstruction of pab, parameterized 
as in Eq. JSJ by the coefficients cy, yielding pab up to 
a local basis rotation [2^ . 

By considering the same four settings at B as 
at A, we can use the linear combination operators 
jA{j) hi Eq. ||TJ and correspondingly js(i) to con- 
struct a complete set of observables, in the basis 
{|l)A|l>fl,|l)A|2) B ,|2)A|l>fl,|2) A |2M, 



3A{j)jB(i) 



(14) 



since the direct products of cr-matrices obey 
tr[(<Tj ®<Ti)(afe®cr ; )] = 4Sj k S u . From Eq. O 
and the relation (ja0)jb(«)) = ^{pAB^j ® &i 
then directly obtain the coefficients Cji as 
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in terms of the measured current correlators and aver- 
aged currents. We emphasize that all elements can be 
determined from sixteen current correlations and eight 
average currents (four at A and four at B). We also note 
that the reconstructed density matrix, due to nonideal 
measurements, might have negative eigenvalues, i.e. it 
might not be positive semidefinite. Schemes to correct 
for this for one and two-qubit states are discussed in e.g. 
ref. H3. 

In the context of two-particle entanglement, it is inter- 
esting to compare the QST-scheme with a Bell Inequality, 
recently discussed for mesoscopic system (see e.g. refs. 
[ill liflj and for a density matrix approach rcf. [HI)- 
Both schemes require the same number of current correla- 
tion measurements. The density matrix reconstructed by 
QST however completely determines the entanglement. 
In contrast, a Bell Inequality can not be used to quantify 
the entanglement |2f|, there are e.g. mixed entangled 
states [2f| that do not lead to a violation of a Bell In- 
equality. 

It is clarifying to illustrate the above scheme with a 
simple example (see Fig. |2J). We consider the Hanbury 
Brown Twiss geometry of ref. [lij . where the number 
of nonzero elements of the one and two-particle density 
matrices are reduced due to the topological properties of 
the conductor. Since no scattering between the upper, 1, 
and lower, 2, leads is physically possible due to the spa- 
tial separation, the one-particle density matrix pa only 
has two nonzero elements, p\\ and p^i- These elements 
are parametrized by cq and C3, obtained by measuring 
(ja(0)) and (m(3)>. 

The two-particle density matrix pab has four nonzero 
elements, pff, p^, P21 an d Pil- Using the relation be- 
tween the coefficients 5,-j resulting from several matrix 
elements being zero, Pab can then be parameterized as 



PAB = 2 
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c 00 + C33 

en + ic 2 i 
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C11 - ic 2 i 

coo - C33 

/ 



(16) 



Conseqently, only four correlations (jA{i)jB(j)) ncc d to 
be measured to completely reconstruct pab , reducing the 
number of actual current correlations needed to be car- 
ried out to 12 for the settings considered here [see Eq. 
1)15(1] . It is interesting to note that in the geometry in 
Fig. |21 considering the tunneling limit for the beam- 
splitters in the source and changing to an electron-hole 
picture |l2j , the emitted state is a true two-particle state 
[l3j . Since the hole currents and current fluctuations are 
directly related to the electron ones, it is possible to em- 
ploy our scheme to reconstruct an electron-hole state as 
well. 

In conclusion, we have presented a scheme for or- 
bital Quantum State Tomography, a complete recon- 



struction of orbital quantum states in mesoscopic con- 
ductors in the transport state. The emitted orbital states 
are manipulated with electronic beamsplitters and de- 
tected with currents and zero frequency current correla- 
tions. With all components experimentally available, our 
scheme opens up for a direct observation of quasiparticle 
quantum states in mesoscopic conductors. 
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